Response maxima in modulated turbulence 
Part II: Numerical simulations 
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Numerical simulations of fully developed turbulence driven by a modulated energy input rate 
or driving force are performed within two dynamical cascade models, the GOY shell model and 
a reduced wave vector set approximation of the Navier-Stokes equation (REWA). The frequency 
behavior of the system response is studied and compared with predictions from a variable range 
mean-field theory, which excludes turbulent fluctuations. In agreement with the mean-field approach 
we find a constant response amplitude for low driving frequencies and a l/w-decay of the amplitude 
for high frequencies. In the mean-field theory, the finite cascade time scale had lead to an oscillating 
behavior of the response amplitude as a function of the driving frequency. In the simulations of both 
models we observe the main maximum. The higher maxima and minima are completely washed out 
by fluctuations. 

PACS numbers: 



I. MOTIVATION 

Many realistic turbulent flows are subject to modu- 
lated driving forces, as e.g. the atmosphere of the earth 
driven by the periodic heating of the sun or the pulsed 
flow through a pipeline. Three dimensional turbulence 
is characterized by an energy cascade from the outer 
length scale, where the forcing acts, to the dissipative 
scale, where most of energy is dissipated, see e.g. 0,0- 
The down-cascading of energy from large to small scales 
takes a characteristic time r. In a statistically stationary 
flow the energy dissipation rate equals the energy input 
rate. In a situation with time dependent energy input, 
on the other hand, this statement will only hold on av- 
erage, whereas the energy dissipation at a certain time t 
is expected to depend on the energy input at an earlier 
time due to the finite time delay of the energy transfer. 

In a previous work Q the effect of an energy input rate 
modulated in time, 



,(t) = e (l + esinut), 



(1) 



with a modulation amplitude e <C 1 and a modula- 
tion frequency u>, has been studied within a variable 
range mean- field theory 0. The response of the system 
can be observed in the second order velocity structure 
function of the flow field at the outer length scale L, 
D L(t) = (((u(x + £,*)- u(x,t)) 2 )) = 6ui, rms , which is 
equivalent to the Reynolds number Re(t) = Ui, rma (t)L/v 
of the flow and the total energy E(t) = ((u 2 )) /2 of the sys- 
tem. Here, u\ >rma is the rms of one velocity component 
and v is the viscosity. The response follows the oscilla- 
tion of the energy input rate with almost constant modu- 
lation amplitude at low frequencies w of the energy input 
rate, whereas the response amplitude strongly decreases 
(oc 1/lo) at higher frequencies. The finite energy transfer 
time r plays a crucial role in this theory. This time r 
is the average time the energy stays within the system 



while it is transported by the interaction cascade from 
the large eddies towards the small eddies, where it is fi- 
nally dissipated. This intrinsic time scale of the system is 
a multiple a of order 1 of the large eddy turnover time T£ , 
corresponding to the sum over the eddy turnover times 
on all scales. r _1 determines the frequency at which the 
crossover takes place between the regime of constant re- 
sponse amplitude and decreasing amplitude. In addition, 
it leads to an oscillating behavior of the system response 
with driving frequency to, where the maxima and min- 
ima are at frequencies connected to the inverse of the 
energy transfer time r. In the limit of large frequencies 
to, the extrema of the response can be estimated to be at 
frequencies u> r ~ n^, n — 1, 2, 3, .... 

Recent experiments on modulated turbulence in a 
cylinder between two counter rotating disks [f| revealed 
evidence for the proposed response maxima. In accor- 
dance with the predictions from the mean-field theory 
for small frequencies a constant response amplitude 
was measured. For large driving frequencies a l/w-decay 
of the velocity response amplitude was observed, again in 
agreement with the prediction from our mean-field ap- 
proach. Note here, that both the velocity response as 
well as the energy response are the same up to a factor 
of 2, in linear order, cf. section VI of reference 0]. 

In the experiments the amplitude of the driving force 
rather than that of the energy input rate is modulated. 
Since the energy input rate is not a controlled quantity 
any more it can serve to measure the response of the sys- 
tem. Of course, also within the mean-field theory we can 
apply a modulated driving force, see 0. The main fea- 
tures, the 1 /w-decay of the energy response amplitude for 
high frequencies and the constant response amplitude for 
low frequencies, pertain. The response maxima are only 
slightly shifted in comparison to the case of a controlled 
and modulated energy input rate. In the case of a mod- 
ulated driving force, as in the experiments, the energy 
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input rate as a response of the system, also shows max- 
ima in addition to the mentioned mean features. These 
are at the same frequencies as the maxima of the total 
energy response amplitude. 

In the mean-field approach, the (intermittent) fluctu- 
ations of the energy and, in particular, of the cascade 
time r are not present. In experiments and numerical 
simulations these fluctuations are of course present, and 
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FIG. 1: Response amplitude A as a function of the driving 
frequency ui for a modulated energy input rate ei n = eo(l + 
esintjt) calculated within the GOY shell model (full dots), 
see section Hi Bl The modulation amplitude is set to e = 0.2, 
and the cascade time delay turned out to be t/tl = a = 
2.54. The stationary Reynolds number is Reo = 7.1 • 10 4 , the 
viscosity v = 1.01875 ■ 10~ 4 , and the large eddy turnover time 
tl = 15.57. Time and length units are set by u, ko and Fo in 
GOY. Our findings are compared with the response amplitude 
as calculated within the mean-field model with the same e and 
r (dashed lines), (a) Log- log plot of the amplitude A versus 
frequency. The long-dashed line denotes the low frequency 
limit of the mean-field theory, A ~ 2/3, and the dotted line 
the high frequency limit, A oc 2/(3o>). The arrow denotes 
ljtl — 1/a — 0.39. Near to this frequency the crossover 
takes place in GOY. Inset: linear scale plot of the response 
amplitude, (b) Log-log plot of the amplitude compensated by 
the asymptotic amplitude, i.e., A/ '(lotl) -1 versus frequency. 
A clear maximum is observed in GOY at a frequency near to 
the maximum of the mean-field amplitude. 



they may lead to broader and less pronounced response 
maxima and minima. Therefore, in this paper we shall 
study the frequency dependence of the response to a mod- 
ulated energy input rate into a system where turbulent 
fluctuations are included. In particular, we shall address 
the question whether the response maxima and minima 
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FIG. 2: Response amplitude A as a function of the driv- 
ing frequency u for a modulated energy input rate ei„ = 
eo(l + esinttrf) calculated within the REWA model (full dots), 
see section Hi! Bl The modulation amplitude is set to e = 0.3, 
and the cascade time scale results to be t/tl = a = 2.94. 
The Kolmogorov constant is found to be b — 83.5 in this 
simulation instead of b exp =6 — 9. The stationary Reynolds 
number is Reo = 1.2 • 10 s , the viscosity v = 5 • 10~ 5 , and the 
large eddy turnover time tl — 0.063. Times are measured 
in units of L^' 6 e^ 1 ^ in REWA. The result is compared with 
the response amplitude as calculated within the mean-field 
model with the same e, r, and b (dashed lines), (a) Log- log 
plot of the amplitude A versus frequency. The dotted line 
is oc The arrow indicates the mean-field crossover fre- 

quency uj^ossTl = (6/b) 3,/2 = 0.019. Inset: linear scale plot 
of the response amplitude, (b) Log-log plot of the compen- 
sated amplitude, i.e., Aj [lotl)~ x versus frequency. A clear 
maximum is observed in REWA at a frequency near to the 
first maximum of the mean-field amplitude. The arrow indi- 
cates the height of the maximum, i.e., a deviation from the 
1/oj-decay by a factor of 1.4 in REWA. 
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can still be well identified in the presence of fluctuations. 
Furthermore, we not only consider a modulated energy 
input rate, but also discuss the slightly different case of 
a modulated driving force in order to compare with the 
above mentioned experiments. 

An appropriate way to numerically study the prob- 
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FIG. 3: Response amplitude A and amplitude of the energy 
input rate A ein as a function of the driving frequency ui for a 
modulated driving force F — Fb(l+esina;t) calculated within 
the GOY shell model (full dots) , see section Ul Gl Fhe modula- 
tion amplitude is set to e/ = 0.2, and the cascade time scale 
is found to be t/tl = a — 2.48. Fhe stationary Reynolds 
number is Reo = 8.6 • 10 4 , the viscosity v = 1.01875 • 10~ 4 , 
and the large eddy turnover time tl = 14.5. Fhe result is 
compared with the response amplitude as calculated within 
the mean- field model with the same e and r (dashed lines), 
(a) Log-log plot of the amplitude A versus frequency. Fhe 
long-dashed line denotes the low frequency limit of the mean- 
field theory, A ~ 1, and the dotted line the high frequency 
limit, A oc 2/(3ojtl). Fhe arrow denotes utl = 1/ffl = 0.40. 
Near to this frequency the crossover takes place in GOY. In- 
set: linear scale plot of the response amplitude, (b) Log-log 
plot of the compensated amplitude, i.e., Aj '(ljtl) -1 versus 
frequency. Fhe dotted line denotes A/ {lutl)^ 1 oc const, (c) 
Linear scale plot of the energy input amplitude A ein versus 
frequency. Fhe mean-field amplitude as well as the GOY am- 
plitude start for low ui with A eiri ~ 1.5 and merge at A ein ~ 1 
for high frequencies. Fhe GOY amplitude shows only the first 
main maximum. 



lem of modulated turbulence would be a direct numeri- 
cal simulation of the Navier-Stokes equation for this spe- 
cific time-dependent energy input rate. However, as we 
need high Reynolds numbers to achieve fully developed, 
isotropic, and homogeneous turbulence and, in addition, 
need the response of the system as a function of time for 
a wide range of driving frequencies, the computational 
demands would be too high. Fherefore, we first study 
the problem within a dynamical cascade model of turbu- 
lence, the Gledzer-Ohkitani-Yamada (GOY) shell model 
[1I1III1EEEE With this model lar S e Reynolds 
numbers and enough statistics within a reasonable com- 
puting time for each driving frequency can be achieved. 
The GOY model has been successfully used in a study 
about decaying and kicked turbulence El- I n addition, 
to be even closer to a numerical Navier-Stokes simula- 
tion and to distinguish between real effects and artifacts 
of the turbulence model, we follow another approach. We 
calculate the response of the system to a modulated en- 
ergy input rate within a reduced wave vector set approx- 
imation (REWA), EE 03, 

where the Navier-Stokes 
equation is solved on a reduced, geometrically scaling 
subset of wave vectors. This method is much closer to 
the Navier-Stokes dynamics than the GOY-model, as it 
contains (i) much more modes than GOY, (ii) it solves 
the Navier-Stokes equation for those modes and not only 
a model equation, and (hi) it is three dimensional. 

Our main results are summarized in Figs. ^ 121 an d 
13 In FigsQ] and [21 the amplitude A of the energy re- 
sponse is shown as a function of the driving frequency for 
both the GOY model (FigHJ and the REWA simulation 
(Fig|2t with a modulated energy input rate. This is com- 
pared with the results of the mean-field model with the 
corresponding parameters, i.e, the same modulation am- 
plitude e and time scale r. In Fig|3|the results from the 
GOY model solutions are shown for a modulated driv- 
ing force and compared with the mean-field model. In 
all cases we observe a constant amplitude for low driving 
frequencies and a l/w-decay for high frequencies. This 
can in particular be observed in the compensated plots 
(parts b of all three figures), where A, compensated by its 
asymptotic amplitude, i.e., Aj '(wr^) -1 , is plotted versus 
frequency. The 1/w-decay of the response means, that 
for fast modulation no response is detectable any more. 
The remaining dissipation rate is that of the stationary 
system itself. 

In the mean-field approach, a sequence of response 
maxima is present for both types of forcing, starting at 
a frequency u) oc 1/r. In the simulations, this main max- 
imum can also be observed, although it is weaker and 
broader, i.e., it is "washed out" by fluctuations. The 
higher order maxima and minima are not visible in the 
simulations, but are completely washed out by fluctua- 
tions. On the other hand, we emphasize that the turbu- 
lent fluctuations in the GOY model are strongly over- 
estimated due to the extreme mode reduction in this 
model. In the REWA simulation, an artificially large 
Kolmogorov constant b is found. Using such large b in the 
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mean-field approach also leads to a considerable weaken- 
ing of the first maximum and a shrinking of the higher or- 
der maxima and minima towards very small amplitudes. 

These results will be explained and discussed in de- 
tail in this paper, which is organized as follows. In the 
next section we study the modulated turbulence within 
the GOY shell model. Before calculating the response of 
the system to a modulated energy input rate as well as 
a modulated forcing in section III Bl and III CI we briefly 
introduce the model and study its stationary properties 
in section Til Al In section ITTT1 we present our findings on 
modulated turbulence within the reduced wave vector set 
approximation. We summarize our results in section llVl 



II. MODULATED TURBULENCE IN THE GOY 
SHELL MODEL 
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A. Stationary properties 

The GOY shell model consists of a set of coupled ODEs 
for one-dimensional complex velocity modes u n |6j, |7J, I& 
EI El El El El These modes u n correspond to velocity 
differences \u(x + r n ) — u(x)\ on scale r n . N modes are 
taken into account, n = 1, 2, N, one complex velocity 
mode per cascade level n, defined by the wave numbers 
k n = X n ko which are equally spaced on a logarithmic 
scale, here, A = 2. The model equations read: 

{it + ^ = U * n+lU * 1 + 2 + 6fc ™-l U T*-i<+l 

-r-cfe n _2tt*_ 1 u*_ 2 ) +F6 n> x, (2) 

where n = 1,...,N, a = 1, b = —1/4, and c = —1/2. 
These are the traditional parameters. We impose bound- 
ary conditions on the u n , i.e., u n = for n < 1 or n > N. 
We use N = 14 shells, a viscosity of v = 1.01875 • 1CT 4 , 
and ko — 2~ 4 . The forcing acts on the largest scale, 
i.e., the first shell, n = 1. F is constant, F = Fo — 
(1 + i) ■ 10~ 2 . Together with v and ko this sets the time 
and length units as well as the Reynolds number. Equa- 
tions (J2J are integrated using a fourth order Runge-Kutta 
scheme with adaptive step size |Tg| . 

With the above chosen parameters the GOY dynamics 
is chaotic 01 • The system is forced on large scales while 
most of the energy is dissipated on small scales. It reaches 
a steady state, in which the velocities are stochastically 
fluctuating. In this sense the system has similar proper- 
ties as three dimensional Navier-Stokes turbulence. The 
scaling behavior of structure functions and dissipation 
has been extensively studied in [13, El O O 123, HH . 
The deviations from K41 scaling due to intermittency ob- 
served in the GOY model are very similar to experimental 
values. Most of the previous studies have been done with 
22 or more shells. As we use here only 14 shells in or- 
der to reduce the computational effort we explicitly check 
some scaling properties of, e.g. the structure functions 
and the energy spectrum in a simulation with constant, 



FIG. 4: Second moment E Ui 2 = {{\Q(u n Un+iu n+ 2 + 
(i^ 1 ) u„-iu n u„+i)\ 2 ' 3 )) (full dots) corresponding to the sec- 
ond order velocity structure function as a function of wave 
number k n = 2 n ko, averaged over a long time interval (t — 
1400tl) and with a stationary, constant forcing, F — Fo- Be- 
tween shell 2 and 9 ISR-scaling behavior, E n ,2 oc k~ a ' 71 (solid 
line) is observed, whereas the shells 12-14 form the VSR. £14,2 
is zero by definition. The dotted line indicates the middle of 
the crossover region between the ISR and the VSR. 



stationary forcing F — Fq in Equations N = 14 then 
turns out to be sufficient. 

The Reynolds number of the system can be defined as 
follows. An outer length scale L is given by the small- 
est wave number k\, L = 1/ki- A typical velocity U 
is the velocity on that scale, ((|ui| 2 ))^ 2 . The average 
((...))t is taken over time. With these length and velocity 
scales the Reynolds number of the present simulation is 
Reo = = 8.6 • 10 4 . The simulated time interval is sev- 
eral hundreds of large eddy turnover times |fciui| _1 . For 
the second order structure function we use the following 
method. In 10] it has been suggested to study the scaling 
of S„ i9 = ((\^s{u n u n+ iu n+2 + u n -iu n u n+ i)\ q / 3 }) in- 
stead of the pure moments of the velocity s n>g = ((\u n \ q )} 
in order to eliminate the period 2 and period 3 oscil- 
lations which are an artifact of the GOY model. The 
second order quantity S nj 2, corresponding to the second 
order structure function, is shown in Fig0]as a function 
of wave number index n = log 2 (A:„/fco). One clearly ob- 
serves an inertial subrange (ISR) between shell 2 and 9, 
where the second order structure function shows scal- 
ing with an exponent near to the K41 value 2/3, i.e., 
S nj 2 oc k~°- 71 corresponding to oc r 71 for the struc- 
ture function. The scaling exponent is not equal to the 
K41 value, because the model shows intermittency cor- 
rections. The higher wave numbers, i.e., smaller scales, 
n= 12 — 14 belong to the viscous subrange (VSR), where 
dissipation takes place. In this range the viscosity term is 
dominant and the velocity decays rapidly with k n , in fact 
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FIG. 5: Energy frequency spectrum of the GOY-system with 
a stationary, constant forcing F — Fo. The spectrum is ob- 
tained by fast Fourier transforming (FFT) an actual time se- 
ries of the velocity u(t) = J2n((-^ e ( Url (^)))) ■ About one fre- 
quency decade of Kolmogorov scaling oc /~ 5 ^ 3 is observed. 
The fc-space the scaling regime of the spectrum is more ex- 
tended, i.e., about three decades, as has been already shown 
in Fig|l 
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FIG. 6: Characteristic time scales of the turbulent energy 
transfer r n (full dots) and of the viscous dissipation t% (full 
triangles) as functions of the level number n. The dashed line 
is a linear fit of log 2 r n , n — 2, 11, and gives r^ 1 * = 22.46 ■ 
2~o.63n rpj^ gnown Tn are obtained by averaging over about 
1400 large eddy turnover times. The dotted line indicates the 
middle of the crossover region between the ISR and the VSR. 



more than exponentially |2fl l2lj . The external forcing of 
the flow acts on shell n = 1, therefore the stirring sub- 
range (SSR) contains the first shell only. The spectrum, 
which is obtained by fast Fourier transforming (FFT) a 
time series of the velocity u(t) = X) ra ((-^ e ( u « W))) an d 
raising it to the power two, can be compared to the en- 
ergy spectrum under the assumption of the Taylor hy- 
pothesis. This spectrum is shown in Fig0for N = 14. 
We observe about one frequency decade of (nearly) Kol- 
mogorov scaling, where the energy decays as oc /~ 5 / 3 
with frequency /. For N — 14 the Kolmogorov scaling 
range is not yet well developed. For less shells this region 
becomes even narrower and the spectrum shows strong 
peaks at some intrinsic frequencies. When reducing the 
number of shells even more, the velocities relax to a sta- 
tionary value without any fluctuations, i.e., the chaotic 
behavior of the system is lost. 

In conclusion, the GOY model with N — 14 shells ex- 
hibits an inertial subrange scaling, although the ISR for 
the frequency spectrum is only narrow, but the spectrum 
in /c-space has a scaling range of about three orders of 
magnitude. This seems acceptable for our goal to study 
the response of the system to a modulated driving, be- 
cause we are only interested in global quantities like the 
total energy E(t) but do not need information on scale 
resolved quantities. 

The time scales in the model have been determined as 
follows. For each shell n an eddy turnover time r„ is 



defined by [Hj: 

This is also considered as the time scale for the turbulent 
energy transfer through the nth level. The time scale 
relevant for the energy loss on level n due to viscosity is 
defined as 

r d = —— (4) 

Both time scales are shown in FigEl 

In the ISR between shell 2 and 9, where the energy 
transfer times are the relevant time scales for the dynam- 
ics, the decrease of the r„ with n is near to r„ cx 2 -2 ™/ 3 
(dashed line in Fig|BJ) as expected for the turnover times 
of eddies of sizes r n /L oc (1/2)™. In this range the dissi- 
pation time scales t% are much larger than the r n , mean- 
ing that the turbulent energy transfer is much faster than 
the viscous dissipation, and therefore the dominant pro- 
cess. In the VSR instead < r n , i.e., on average the 
energy is dissipated by viscosity before it can be trans- 
ferred to smaller scales. 

The largest eddy's turnover time is in general de- 
fined by the velocity on the outer length scale L, i.e., 
on the length scale of the forcing, which in this case 
is 1/fci. However, in this model, the time t\ is dis- 
turbed due to finite size effects. Therefore, we extrap- 
olate from the turnover times of the other shells. A lin- 
ear fit in FigEl for log 2 r n with n — 2,..., 11 leads to 
r[ lt = 22.46 • 2~ 63n . Without intermittency one would 
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have t„ oc 2~ 2 ™/ 3 ; the small deviation corresponds to the 
intermittent scaling of u n oc fc~ ' 37 or an intermittency 
correction of <5£i = £1 — 1/3 ~ 0.04. The extrapola- 
tion for n — 1 yields for the large eddy turnover time 
t l = r( lt = 14.5. The time scale corresponding to the 
energy transfer time r used in the mean-field model 0] is 
the sum over the eddy turnover times of all energy-input 
and inertial-range shells, i.e., here r ~ En=i T ™ — 35.9. 
The factor t/tl = a between the transfer time and the 
large eddy turnover time is then a = 2.48. 



Constant energy input rate 

Until now, we have considered a constant forcing F = 
Fq. The resulting energy input rate e, n (i) = ((u*(t)Fa)) 
then fluctuates around its mean value because of the 

(^-fluctuations. In the mean- field theory Q| the en- 
ergy input rate ei„ is constant instead. For closer com- 
parison we also consider another type of forcing in the 
GOY-model: F(t) = e Ui{t) /\m{t)\ 2 . This forcing F(t) 
fluctuates as ui(t). Then the energy input rate is ei„ = 
((u*F)) = £o = const by definition. The ISR-scaling be- 
havior as well as the energy spectrum then turn out to be 
similar to the previously discussed ones with the constant 
forcing F = F . The energy transfer time r is slightly 
larger in this case, namely t = 39.5, and the large eddy 
turnover time is tt, — r[ lt — 15.57. Again, the large eddy 
turnover time is extrapolated from t 2 , t 1:l . This leads 
to the factor t/tt, = a = 2.54 between the total time 
delay of the energy cascade and the large eddy turnover 
time. In the following sections we will study the time- 
dependent cases where either the energy input rate ej n , 
i.e., F = €qUi/\ui\ 2 ( section Til Bll or the forcing F = F 
is modulated (section III Cf) . 



B. Modulated energy input rate 

In this section we apply a modulated energy input rate 
to the GOY model, i.e., we set the forcing F = F(t) in 
equations (J2J as 

Til 

F(t) = e 0l —^(l + e sin wt) (5) 

with a modulation amplitude e = 0.2. Then, the result- 
ing energy input rate ej„ is 

e m (t) = ((ul(t)F(t))), (6) 
= eo(l + e sinwi), 

and has a prescribed modulation amplitude eeo by defi- 
nition. The total energy of the system 

^) = ^EK(*K(i)». (7) 

n=l 




10 20 30 40 10 20 



t/x L t/t L 

FIG. 7: Energy input rate e£? rm = e in /eo (dotted lines) and 
energy content E norm (solid lines) for four different modula- 
tion frequencies ui calculated in the GOY model. The energy 
input rate ei n is modulated with a modulation amplitude of 
20% of the constant energy input rate eo, e = 0.2, according 
to Eq.JSJ. Also included is the fit according to Eq.@ for the 
energy E norm (dashed lines, indistinguishable from the solid 
lines), (a) lut l = 0.0151, (b) ujt l = 0.151, (c) ujt l = 0.787, 
(d) ujtl = 2.025. 

is calculated for a wide range of driving frequencies u> 
in order to study the frequency behavior of the response. 
The brackets ((...)) denote the ensemble average. This en- 
semble average is performed as follows. From a long sta- 
tionary simulation we collect an ensemble of 1500 start- 
ing configurations which we then let evolve according to 
Eqs.J2J but now including the modulation of the forcing 
F(t), Eq.JSJ), and average over these 1500 time series. To 
ensure that the different realizations can be considered 
as statistically independent, the time delay between the 
successive starting configurations for the different real- 
izations is chosen to be about 100 large eddy turnover 
times. The adaptive step size routine controlling the 
Runge-Kutta integration algorithm for ODEs does not 
produce the same time steps for all time series. To over- 
come this, we have calculated also the same number of 
equidistant time steps by spline interpolation for all time 
series. 

The oscillating response of the system A(£) is then 
studied in terms of the ratio between the energy E(t) 
with modulated energy input and the energy E$ (t) with- 
out modulation, namely 

E norm (t) = ^ = l + A(t). (8) 

E and Eq both are averaged over 1500 realizations. In 
spite of the averaging not only E but also E a still contains 
(weak) fluctuations. Therefore, we write E (t), as E is 
still slightly fluctuating around its mean value. E and Eq 
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contain about the same size of fluctuations. Accordingly, 
the energy input rate ej n is normalized by its stationary 
value, e?° rm (t) = e in (t)/e = 1 + esinut. 

In Fig[7| the input rate e™° rm (t) and the energy 
E norm (t) are plotted for four different driving frequen- 
cies. For the two low frequencies where ujtl <C tl/t = 
1/a ~ 0.39, the energy follows the oscillation of the en- 
ergy input rate with almost constant, but smaller ampli- 
tude. For higher frequencies the amplitude of the devia- 
tions of the normalized energy from its stationary value 
1 strongly decrease, and a phase shift with respect to the 
energy input becomes visible. The same behavior of the 
energy has been observed in the mean- field theory |3j. 

To quantitatively access the frequency behavior of the 
response amplitude, we calculated time series of the total 
energy E(t) for 85 different driving frequencies varying 
over almost 3 decades between 0.012 < lotl < 9.3. The 
chosen frequencies are approximately equally spaced on 
a logarithmic scale. The normalized energy E norm (t) is 
fitted by a function of the form 



E norm (t) = E const + eAsm {cut + $), 



(9) 



with three free parameters: E const , the amplitude A, and 
the phase shift $. E const is near to 1 for all frequencies, 
i.e., E const = 1.0022 ±0.0032. The fits © are included in 
FigQas dashed lines but they are mostly indistinguish- 
able from the solid lines for the energy itself. 

Figd°f section[I]shows the amplitude A, resulting from 
the fitting procedure, as a function of the dimensionless 
frequency lut^. A is almost constant for low frequencies 
and has a value of about 2/3. For higher frequencies the 
amplitude decreases as oc 1/ui. The same features have 
been observed in the mean- field calculations, see dashed 
lines. The long-dashed line in Fig^ represents the low 
frequency limit of the mean- field theory, A ~ 2/3, and 
the dotted line the high frequency limit, A oc 2/(3wtl). 

The crossover between the regime of constant ampli- 
tude and the one of l/w-decay of the energy response 
takes place at t-u cross TL — 1/a ~ 0.39, which is indicated 
by the arrow in Fig^ In the mean-field approach this 
crossover is always at u>tl — 1, independent of the fac- 
tor a between the large eddy turnover time and the total 
time scale of the energy transfer. In experiments 5] the 
crossover frequency has been used to measure the energy 
cascade time scale. The present simulations confirm that 
this frequency gives the correct order of magnitude for 
the energy transfer time. 

Response maxima, as observed in the mean-field model 
at frequencies connected with the inverse energy trans- 
fer time, are difficult to be identified in Fig^,. There 
is some structure visible at ojtl ~ 0.31 and lotl — 1.57. 
In Fig^D, where the amplitude A, compensated by the 
asymptotic amplitude (lotl) , i.e., A/ (ujtl)^ 1 is plot- 
ted versus frequency, this structure becomes more ev- 
ident, and we see a clear maximum at a frequency of 
about ljtl — 1.57. This maximum probably corresponds 
to the mean-field maximum. Of course, the maximum in 
GOY is broadened and weakened due to the large fluc- 



100 



Q. 




20 



100 



200 



t/T L 



FIG. 8: Time series of the cascade time r(t) for non- 
modulated forcing F = eo^i/l^il 2 in the GOY model. Strong 
fluctuations are observed. Inset: Probability distribution of 
t/t l . The mean is a = 2.54, and the width FWHM/2=0.46a 
giving about 50% fluctuations. 



tuations, and the higher order maxima and minima are 
apparently washed out completely. As in the mean-field 
theory no fluctuations are included, the energy cascade 
time r is considered to be constant. However, in the 
GOY model this assumption is not true, as can clearly 
be seen in Fig|8] Here, a time series of the cascade time 
r(t) = 53 n =i T n(t) is plotted, computed within the GOY- 
shell model with non-modulated forcing F = eoUi/|iti| 2 . 
The inset shows the probability distribution of t/tl,. 
This distribution has its maximum at t/tl = 2.39, al- 
most at the mean cascade time a — 2.54, and a width 
(FWHM/2) of about 0.46a. The width is almost half 
the size of the mean which indicates that the transfer 
time fluctuates strongly and therefore we have to expect 
that the response maxima are more or less washed out. 
However, these strong fluctuations are considered as an 
artifact of the GOY-model and not as a feature of real 
turbulence. The GOY-model contains only one velocity 
mode per cascade level instead of infinitely many modes 
in real turbulence. This one-mode approximation leads 
to an overestimation of the fluctuation strength. In order 
to confirm this, we performed another simulation with 
more modes per level within the reduced wave vector set 
approximation (REWA) of the Navier-Stokes equation. 
This will be presented in section ITTll 



C. Modulated driving force 

In this section we present further results within the 
GOY model based on a non-fluctuating driving force Fq 
which is regularly modulated as was eo in the previous 
section. This case may be more comparable to the exper- 
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imental method in reference [5j , because there the driving 
force is modulated. What cannot be modeled with GOY 
is the spatial inhomogeneity in the experiments 
In equations J2J) we now apply a forcing 

F = F(t) =F (l + e/sinwt), (10) 

with a modulation amplitude of e/ = 0.2. As in section 
III Bl we calculate the ensemble averaged time series of the 
energy input rate ei n (t) , see first line of equation © , and 
the total energy of the system E{t) cf. equation Q, for 
89 different frequencies between 0.0144 < ljtl < 3.04, 
again logarithmically equally distributed. The normal- 
ized energy E norrn (t) and energy input rate ef° rm (t) are 
then fitted by a function according to equation with 
the parameters E const , A, and e in ^ consU A ein , $ eiri , 
respectively. 

The amplitudes A and A eiri are plotted in Fig|3] of 
section [I] as a function of the dimcnsionlcss frequency 
lotl. Also in this case of a modulated force, the re- 
sponse amplitude is almost constant for small frequen- 
cies, namely A ~ 1, and decreases as 1/w for high fre- 
quencies, see FigOt (full dots). Again, the long-dashed 
line represents the low frequency limit of the mean-field 
theory for A (which is A ~ 1 in this case) as well as 
the dotted line the high frequency limit. As in the pre- 
vious section, the crossover frequency between the two 
regimes is determined by the energy transfer time, i.e., 
Ucrossi~L — 1/a ~ 0.40 (as a = 2.48 in this case), which is 
marked by the small arrow in FigO^- The amplitude of 
the energy input rate A ein starts with a value of about 
1.5 for low frequencies and merges towards 1 for high fre- 
quencies, see FigEfc- This indicates that at very large 
frequencies the velocity is not oscillating any more as it 
only feels a mean constant force. The oscillations of the 
energy input rate are then only a consequence of the oscil- 
lation of the driving force F. In the mean-field theory we 
have observed the same trend for both amplitudes. The 
corresponding mean-field results are included as dashed 
lines. 

Both amplitudes A and A ein show a maximum at a 
frequency near to the crossover frequency ujtl — tl/t = 
a^ 1 = 0.40. In the compensated plot FigO), where 
A/^tl)^ 1 is plotted as a function of frequency a clear 
deviation from the dotted line representing A oc I /to can 
be observed. At this frequency, the mean-field theory 
predicts a first maximum for the energy response am- 
plitude (Fig|3K) and a maximum directly followed by a 
minimum for the energy input rate, see Fig|3|;. This fre- 
quency is connected to the energy transfer time r. In the 
mean-field model further maxima of the energy response 
and wiggles of the energy input rate are observed at mul- 
tiples of this frequency. However, in the GOY model all 
further maxima and minima are washed out presumably 
because of the strong fluctuations. The fluctuations of 
the energy cascade time scale are found to be similar as 
in the case of a fluctuating force shown in FigJH] in the 
previous section. 



III. MODULATED TURBULENCE IN A 
REDUCED WAVE VECTOR SET 
APPROXIMATION 

As was pointed out in the introduction of this chapter 
a full numerical simulation of the Navier-Stokes equa- 
tion for modulated turbulence at high Reynolds numbers 
is still not possible or requires low Reynolds numbers. 
Therefore, we first have considered the GOY shell model. 
This model correctly describes many features of turbu- 
lence, however, due to the one-mode approximation in 
each cascade level, it contains various artifacts. Namely, 
it strongly overestimates the strength of the fluctuations. 
The aim of this section is to study the problem of modu- 
lated turbulence within another model, the reduced wave 
vector set approximation (REWA) 0, 0, 0, which 
is much closer to the Navier-Stokes equation than the 
GOY model and contains much more modes per cascade 
level. Of course, as compared to full numerical simu- 
lations of the Navier-Stokes equation, it still contains a 
mode reduction in order to make the computational ef- 
fort reasonable for the desired high Reynolds numbers. 
The present approximation has been introduced and ex- 
tensively studied in 0, 0, 0] . Here, we use it to- 
gether with a time-dependent driving. For completeness 
we briefly explain the approximation before we present 
the results with modulated driving. 

A. The reduced wave vector set approximation 

The velocity field u(x,t) is Fourier transformed into 
plane waves, u(x,t) = ^]pi((p, ije'P' 1 . Periodic 
boundary conditions are applied on a periodicity vol- 
ume (2ttLq) 3 . The wave vectors p are given by p — 
(Pi) = (^Lq 1 ), with rii = 0, ±1,±2, .... In order 
to efficiently deal with the large number of modes in- 
volved, the reduced wave vector set approximation se- 
lects a limited number of modes by admitting only a 
geometrically scaling subset K = (J ; Ki of wave vec- 
tors, i.e., u(x,t) = ~}2p eK u(p,t)e l P' x . On this subset 

K = {pn\n = 1, N, 1 = 0, l m ax} the Navier-Stokes 
equation for incompressible flow, 

= - V ( p ii))* u .( p v)) + fi ( p m) (ii) 

-i M ijk (p«) Yl «j(«i)«*(9a). 

q lt q 2 eK,q-,+q 2 =p^ 

together with the continuity equation, 

pW.„(pW) = , (12) 

is solved. My/j is the coupling matrix, My&(p) = 
\(Pi P ik(p) +PkP l r j(p)), where (p) is the orthogonal 
projector to p. The subset K consists of a basic subset 
A'o = {pi , n = 1, N} together with its scaled replicas 
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pg>= 2^,1 = I,. 



■ ilmax- In the present simulation 



we take N = 74 wave vectors p 



(0) 
n , 



K = {±(2, 2, 2), ±(-1,2, 2), ±(-2, 1,1), ±(3, 0,0), 
±(4, 1, 1), ±(4, -2, 1), ±(-3, 3, 3), ±(-5, 1, 1), 
±(4, 4, 1), ±(3, 3, 0), ±(1, 1, 10), ±(-10, 5, 5) 
+permutations} (13) 

These wave vectors {J t Ki are chosen such that they span 
a wide range of length scales, but still dynamically in- 
teract to a good degree. For the Navier-Stokes equa- 
tion this means that as many Navier-Stokes interactions 
p = q 1 +q 2 as possible between the wave vectors in K are 
allowed for. The largest eddies of the order Lq are rep- 
resented by the wave vectors in Ko, whereas the subsets 
Ki contain wave vectors of smaller and smaller eddies. 
The choice of the smallest eddies, i.e., the value of l ma x, 
depends on the kinematic viscosity v. l max and v are 
adjusted such that the velocity amplitudes u(pn ma!C \t) 
of the smallest eddies are almost zero. In this simulation, 
v is chosen as^ = 5-10~ 5 and the number of levels as 

lmax ±1 9. 

To maintain the turbulent flow we apply a forcing as 
in|lil3: 



f( p ,t) = l e °E^S™^F( 1 + esin ^' peKl ' 
\ 0, P^K„ 



(14) 



The subset Ki n of Kq by choice contains the 
wave vectors with the three smallest lengths. 
Ki n contains 14 vectors, namely Ki n = 
{±(2, 2, 2), ±(-1, 2, 2), ±(-2, 1, 1), +permutations}. 
In reference [l6|| it has been shown that the statistics of 
the solutions of the equations of motion do not depend 
on the particular choice of Kq. This forcing corresponds 
to the same type of forcing, which has already been 
applied to the GOY model in section IIIBI It enforces 
the energy input rate to be modulated: 

e in = (( J2 u *(Pn) ■ /(Pi ))) = + esin ^)' ( 15 ) 
p^eK 

Eq.JTll) is a set of 3N(l ma x + 1) coupled ODEs for the 
complex mode amplitudes Wi(pi^) which is numerically 
solved within the Burlisch-Stoer integration scheme with 
adaptive step size [l8j. Length scales are measured in 

units of Lq and time scales in units of Lg^ 3 e 1 ^ 3 . A 
Reynolds number can be defined as follows: The wave 
length A of the smallest wave vector gives an external 
length scale L — 2n/^/6, and a typical velocity on that 
scale is determined by the rms of one velocity component, 
ui.rms- Then, in our case, the Reynolds number is Re = 
Ul - rmsL — 1.234 ■ 10 5 because from the simulations we 
obtain ui. rms — 2.405. 

The main features of fully developed turbulence as ir- 
regular velocity signals, characteristic scaling of structure 



functions, etc. are well described within this approxima- 
tion, as has been shown in EEEEl. The REWA solu- 
tions show small scale intermittency, which is produced 
by the competition between down-scale en ergy tr ansport 
and viscous dissipation on the small scales |16ll22| . Other 
mechanisms leading to intermittency in turbulence as e.g. 
nonlocal interactions between wave vectors are underes- 
timated in this approximation |23| . The down-scale en- 
ergy transport in the REWA fluid is less effective than in 
real turbulence, because in this approximation the larger 
wave vectors are more and more thinned out [24j. This 
is in contrast to the case of the complete set of wave vec- 
tors (e.g. in full grid simulations) where the density of 
states increases oc p 2 , whereas in the reduced wave vec- 
tor set K the number of admitted wave vectors decreases 
as 1/p 01 ■ I n reference [24| it has been shown that 
this reduced energy transport leads to an overestimation 
of the Taylor Reynolds number of the system as well as 
the Kolmogorov constant b, defined by D(r) = 6(er) 2 / 3 , 
by roughly one order of magnitude for our choice of N. 
In the present simulation we obtain b = 83.5 instead of 
b = 6 - 9 as in experiments [M EE El Since D(L) 
is the energy density oc ((u 2 )) of the fluctuations in the 
fluid system, the large b value indicates that even in the 
REWA approximation the strength of the fluctuations is 
highly overestimated. The large Kolmogorov constant 
will change the relevant time scales in the system, as will 
be shown in section IlII Bl 

The characteristic time scale for the turbulent energy 
transfer on scale I can be estimated as: 



r(l) = 



P Urm 



(16) 



where denotes the mean wave number on scale Z, i.e., 
it is the mean inverse eddy size in K\. As in the GOY 
model, Eq.(0}, the time scale of viscous dissipation is 
T d = [^(p^) 2 ] -1 • Again, in the ISR t{1) > t^, whereas 
in the VSR > t(Z). From a simulation with stationary 
forcing, i.e., e = in Eq. (|T4*|) . the time delay of the 
energy down-transport r sum is then estimated by the sum 
of all t{1) in the ISR, r sum = J2ieiSR T ( l ) - °- 186 - The 
largest of these t(1), on the largest scale, t(0) = 0.0632 
can be regarded as a large eddy turnover time t/,. Thus, 
T S um — 2.94t£, and the factor between the cascade time 
scale and the large eddy turnover time is T stim /Tx = a = 
2.94. 

As we have seen in the GOY model, FiglSI the energy 
transfer time is strongly fluctuating. We attributed these 
strong fluctuations to the one-mode per level approxima- 
tion of the GOY model. Fig|5] shows a time series of the 
energy transfer time T sum (t) (use Ea. (|16fl with vS l \i)) 
in the present reduced wave vector set approximation 
(REWA), and, in the inset, the distribution of this time 
scale. Clearly, the fluctuations are much weaker than in 
the GOY model; they are zero in the mean-field approxi- 
mation. The distribution is centered around T sum /rL = a 
with a width (FWHM/2) of about 0.02a. 
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FIG. 9: Time series of the cascade time T sttm (i) in REWA. 
Inset: Probability distribution of T sum /TL- The fluctuations 
are considerably smaller than in the GOY model. Note 
the different scales in this figure and Fig|8] The mean is 
Tsum/r L = a = 2.94, and the width FWHM/2 = 0.02a giv- 
ing about 2% fluctuations. 



B. Modulated energy input rate 

The response of the system to a modulated driving 
force, cf . Eq. ifTil) , is calculated now in terms of the total 
energy of the system 



z=z max 

£(*) = 2« E E i«(p.')i 2 »- 



(17) 



The modulation amplitude of the energy input rate ej n 
fEa. (|T3|l % ) is chosen as e = 0.3. The average ((...)) is per- 
formed as follows. We average over 25 to 30 realizations, 
which are obtained from Eas. (|ll|) with different starting 
values. The delay between the different starting values is 
one driving period. For the higher frequencies the period 
lixjui becomes too small to ensure statistical indepen- 
dence of the different realizations. Then, we chose the de- 
lay between the successive starting values of the different 
realizations to be several driving periods such that it is at 
least 60r stlJ „. The response is calculated for 150 (approxi- 
mately equally spaced on a logarithmic scale) frequencies 
between 0.00016 < ojt l < 3.0. The energy is normalized 
by -Ecu calculated from a stationarily forced solution with 
e = and averaged as E. The oscillating response A(t) 
of the system is defined in the same way as for the GOY 
model, see Eq. (JHJl . Then, as for the GOY calculations, 
the averaged and normalized signals E norrn (t) are fitted 
with Eo-.J^J. The fit parameter E const is again near to 1 
for all frequencies, E const = 1.0064 ± 0.0065. In FiglTUI 
the time averaged responses and the normalized energy 
input rates are plotted for four different driving frequen- 
cies. Also the fits according to Eq.® have been included 



FIG. 10: Energy input rate egf™ = 1 + esinut (dotted lines) 
and energy E noTm (solid lines) for four different modulation 
frequencies ui as calculated in the REWA simulation. The 
energy input rate is modulated with a modulation amplitude 
of 30% of the constant energy input rate, i.e., e = 0.3 in 
Eg. (1141 . Also included is the fit to the energy data cf. Eq.© 
as dashed lines but these are in all cases indistinguishable 
from the solid lines. The averaged time series of E norm are 
repeated once for better visibility, (a) ujtl = 3.16 • 10 -4 , (b) 
ujt l = 3.16- 10" 3 , (c) ujtl = 3.16- 10 -2 , (d) ujtl = 0.316. For 
larger ujtl the energy is indistinguishable from 1 on this scale. 
The crossover to the l/o;-decay regime is in this simulation 
at ujtl — 1-1 • 10 -2 , between the frequencies of b) and c). 



as dashed lines but are indistinguishable from the solid 
lines for the energy signal itself. 

We observe in Fig^|for REWA the same features as in 
FigQfor the GOY model and in Fig.l of reference Q for 
the mean-field model. For the two lower frequencies the 
response amplitude remains almost constant and is about 
2/3 of the amplitude of the energy input rate, whereas 
for the two higher frequencies the response amplitude 
strongly decreases. This trend becomes more clear in 
FigEk m section where the amplitude of the response 
- determined from the fit J^J) - is shown as a function of 
the driving frequency (full dots) . For low driving frequen- 
cies A ~ 2/3, whereas for high frequencies the amplitude 
decreases as 1/ui. The crossover between the regime of 
constant amplitude and the one of 1 /w-decay takes place 

at UJ r 



S T£ ~ 0.011, i.e., at a much smaller frequency 
than expected from the original case of the mean-field 
theory in which the crossover was at oj^.^ ss t l = 1 with 
b — 6. We understand this as follows. In section till Al 
it was mentioned that the Kolmogorov constant in the 
REWA simulation is b = 83.5 instead of b = 6 — 9 as in 
experiments. In the figures of the mean-field approach Q 
we have set 6 = 6. The mean-field solution for a general 
b revealed that the crossover frequency decreases with in- 
creasing b while the positions of the response maxima are 
left unchanged. For b — 83.5 the mean-field crossover fre- 
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quency is at ^>cross T L — (6/b) 3 / 2 ~ 0.019 in close agree- 
ment to what we observe in the REWA simulations. The 
response amplitude calculated from the mean field model 
with b = 83.5 is included in Fig[21} as dashed line. Apart 
from the changed crossover frequency we observe that, in 
the mean-field calculations, the first response maximum 
at lotl ~ 0.1 is considerably smaller and broadened as 
compared to the case with 6 = 6, given as the dashed 
line in Fig^,. In agreement with this, our REWA simu- 
lations (with a value b = 83.5) show a broad maximum 
in the response amplitude at ujtl ~ 0.028. This means 
it occurs at a similar frequency as the mean field model. 
The maximum becomes more clear in the compensated 
plot, FiglSb, where A/ (ujtl)^ 1 is shown as a function of 
frequency. There, we observe a deviation from the 1/oj- 
decay of the amplitude by a factor 1.4 at the maximum 
in the REWA simulations. The mean field maximum has 
a height of 2.8. The subsequent maxima and minima 
in the mean field model occur at frequencies where the 
amplitude is already very small {A < 10~ 2 ) because the 
crossover to the 1 /w-regime takes place at a much smaller 
frequency whereas the response maxima stay at the same 
frequencies as for a smaller b. Therefore, the higher order 
maxima are not visible in the REWA simulations. The 
cascade time shows about 2% fluctuations in REWA as 
shown in Fig. [§] However, at small response amplitudes 
these fluctuations are already large enough to wash out 
the higher order response maxima. 

In conclusion, the REWA system reproduces qualita- 
tively the features of modulated turbulence as predicted 
by the mean-field model including the first response max- 
imum. The latter is considerably weakened due to the 
large Kolmogorov constant in REWA. Another conse- 
quence of the large b is, that the crossover between con- 
stant amplitude and 1 /w-decay is shifted towards smaller 
frequencies, and therefore the higher order maxima and 
minima are already at very small amplitudes where the 
fluctuations in the cascade time scale are finally large 
enough to wash them out. We cannot clarify, at present, 
how close the response in direct numerical simulations 
- which lead to an order of magnitude smaller b and 
thus have much smaller fluctuations - will come to the 
mean-field features, but we expect a clearly visible first 
maximum at least. 



IV. CONCLUSIONS 

We have simulated the response of modulated turbu- 
lence within two numerical models. Namely, we have 
used the GOY shell model and the reduced wave vec- 
tor set approximation of the Navier-Stokes equation 
(REWA) . The results are compared with predictions from 
a mean-field theory. For a modulated energy input rate 
this mean-field theory had predicted a constant response 
amplitude for low frequencies and a l/w-decay for high 
driving frequencies. In addition, at certain frequencies 
connected with the energy cascade time scale, a sequence 



of maxima and minima of the response amplitude is ob- 
served. 

Both numerical models well reproduce the basic trend, 
i.e., the constant amplitude for small uj and the 1/uj- 
decrease for large to. The main response maximum can 
be observed in both numerical models, although it is 
weakened due to fluctuations. The higher order max- 
ima and minima as predicted by the mean-field theory 
cannot be identified in the simulations. Obviously, they 
are strongly washed out by fluctuations. In the GOY 
model the large fluctuations are explicitly visible in a 
broad cascade time distribution. We attribute these fluc- 
tuations to the strong mode reduction in the model, i.e., 
they are an artifact of the model properties and not a 
feature of real turbulence. In the reduced wave vector 
set approximation of the Navier-Stokes equation these 
fluctuations are much weaker, and, we believe, more re- 
alistic for real turbulence, as more modes are taken into 
account. However, due to an overestimated Kolmogorov 
constant b in REWA the higher order maxima are consid- 
erably reduced, and therefore washed out by the fluctua- 
tions although being smaller than the GOY-fluctuations. 
Therefore, we believe, that in real turbulence with a real- 
istic Kolmogorov constant the first maximum should be 
clearly observable and possibly also the higher order ex- 
trema in the response. Thus, the predictions of the mean- 
field model, which excludes all fluctuations, might be 
quite reasonable for real turbulence. To further study the 
response maxima numerically, it is necessary to perform 
full numerical simulations of the Navier-Stokes equation, 
as then all relevant time scales including their fluctua- 
tions are reproduced realistically, which turned out to be 
essential for the observation of the response maxima. 

Recent experiments on modulated turbulence [SJ re- 
vealed evidence for the response maxima. These experi- 
ments may be more comparable to the case of a modu- 
lated force instead of a modulated energy input rate. We 
have studied this case also within the mean-field model 
and have found basically the same behavior of the en- 
ergy response as for a modulated energy input rate. In 
addition, the amplitude of the energy input rate showed 
"wiggles" at the same frequencies where the energy re- 
sponse had maxima. In the experiments the response 
maxima were measured in the energy input rate, which 
can be regarded as a response of the system as well in 
this case. Also the constant amplitude for low driving 
frequencies and the 1 /w-decay of the velocity response - 
which in leading order is corresponding to a 1 /w-decay of 
the energy response, as well - have been observed in the 
experiments. Here, we have studied the case of a modu- 
lated driving force within the GOY shell model. Also in 
this simulation, the response amplitude behaves basically 
as in the mean-field model, i.e., it decreases as l/u>. The 
energy response amplitude as well as the amplitude of 
the energy input rate show the main maximum. Due to 
the above mentioned large fluctuations all higher order 
maxima are washed out in the GOY model. 

There are two regimes in the frequency behavior of 
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the response amplitude, namely a constant amplitude at 
low uj and a decreasing amplitude at high uj. The present 
simulations give further confidence that the crossover fre- 
quency between these two regimes gives the correct or- 
der of magnitude of the cascade time scale, i.e., in ex- 
periments it can be used to measure this time scale as 
suggested in reference 

Both models in the present study were able to repro- 
duce the main features of the frequency behavior of the 
response amplitude in modulated turbulence as predicted 
by the mean-field model, however, both also have their 



shortcomings, which prevent us from correctly predicting 
the behavior of real turbulence in all quantitative details. 
Therefore, we believe, that it is worth to further study 
modulated turbulence numerically as well as experimen- 
tally. 
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